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ABSTRACT 

The development of a formalism for the magnetic and thermal prop- 
erties of the Heisenberg ferromagnet in an applied magnetic field is pre- 
sented. By the use of Green's function theory the magnetization and 
susceptibility are determined in the random phase approximation. It is 
shown that the thermodynamic quantities can be expressed in analytical 
form for all temperatures, both above and below the Curie point. 

The object of this work is to evaluate the use of ferromagnetism in 
a refrigeration system. From the analytical results on magnetization 
and susceptibility, equations describing adiabatic and isothermal proc- 
esses are presented. Magnetization and entropy are presented as func- 
tions of temperature and field. 


Keywords: Green’s function. Thermodynamics, Eptropy, Refrigeration, 
Critical point, Magnetocaloric effect 
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INTRODUCTION 

According to the Weiss theory a ferromagnetic body possess an in- 
ternal field which is proportional to the magnetization. As the tempera- 
ture of a ferromagnet is increased, the magnetization decreases until, 
at a temperature known as the Curie temperature T^, the sample be- 
comes paramagnetic. The process by which this change from ferro- 
m^netism to paramagnetism occurs is referred to as the ferromag- 
netic phase transition. Investigations have shown that electrical, me- 
chanical, and many thermodynamic properties of a material are altered 
when the material undergoes a phase transition. 

The effect of an external magnetic field is twofold: (1) the mag- 
netization is increased somewhat above its zero field value and (2) the 
critical point ceases to exist. Even though temperature tends to de- 
stroy this alinement, the field causes some ordering to be present. 

Thus instead of an abrupt disorder at the critical temperature there is 
a more gradual transition. 

This report is concerned with the investigation of the thermodynamic 
properties of the Heisenberg ferromagnet in an external magnetic field. 
The quantities studied are the magnetization, susceptibility and the de- 
rivative of magnetization with respect to temperature. The last of these 
determines thermodynamic quantities, such as entropy and adiabatic 

demagnetization (magneto - caloric effect), which are useful for inves- 

1-2 

tigating various recently proposed refrigeration systems. 

MODEL 

The Heisenberg model is based on a solid where the magnetic elec- 
trons are in states localized about the lattice sites with exchange inter- 
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actions between electrons taking place between nearest neighbor pairs. 

The model does not take into account an itinerant electron picture. The 
model is thought to be very good for insulating ferro magnets such as 
EuO and poor for conductors such as iron and nickel. However, it ap- 
pears that this model gives better results for conductors at low tem- 
peratures than those calculated from a band theory approach as shown 
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by Arglye, Char ap and Pugh. It should be noted however, that the 
Heisenberg model with its surprising success does not take into account 
the spreading of the electronic energy levels into bands. 

The Heisenberg ferro magnet with spin 1/2 was analyzed by Bogolyubov 
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and Tyablikov using the techniques of double-time, temperature depend- 
ent Green’s functions. A convenient review of Green’s functions and 

Tyablikov’s application of them to ferromagnetism is given by D. N. 
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Zubarev and references contained therein. 

The Hamiltonian for the Heisenberg model is 


H = - 2 Ji, ifesf . 1 (s*s: . S-S-)] 

i i,J ^ 


( 1 ) 


where 

Sf = Sf ± iSj (2) 

J J >1 

is the Bohr magneton, g the Lande g factor, H^ the applied magnetic 
field which is assumed to be along the z-direction, S. is the spin oper- 
ator for a spin at site i, and where J. . is the exchange interaction 
between spins on sites i and 3 . The sum is carried over all sites in 
the crystal. The exchange interaction is assumed to be a function only 
of the distance between sites. The self exchange terms such as J.. or 
J.. vanish. 

J J 
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CALCULATION OF THERMODYNAMIC QUANTITIES 

Since the thermod 3 miamlc quantities of interest are the magnetization 

and entropy of the system, one is interested in evaluating the correlation 

functions of the form /s"S'L ) and hence of the Green’s function (/s^; S' )) 

\ g m/ \\ g’ mil 


one can derive the 


starting with the Green’s function ®m)) 
correlation function in the random phase approximation. The results 
are 

-1 




- 1 


(3) 


where 

Ek=2{s‘^)j(0) 

and 


■ J(k) 
J(0) 


+ figgHo 


(4a) 


m 

The sum in equation (3) goes over all N lattice vectors in the first 
Brillouin zone. 

Application of a magnetic field is assumed sufficient to orient 
the net magnetization along the direction of the field. The field also in- 
troduces additional long-range order. The long-range order is thus due 
to the internal field and to the applied magnetic field. The average z 
component of the spin is a measure of this long-range order. The ther- 
mal average ^S^^ is proportional to the magnetization. 
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To obtain the magnetization one calculates by using the relation 
(S^) = — - (5) 

^ ' 2[1 + 2cp] 

where 



The sum in equation (6) must be evaluated over all values of k in the 


first Brillouin zone of the appropriate lattice. Except at the very low- 
and high-temperature limits numeric^ methods are usually used. Such 
numerical solutions are, however, somewhat difficult. One of the pur- 
poses of this paper is to show that analytical solutions of equation (6) are 
possible. 

Consider only crystals with cubic symmetry such as bcc. One can 

replace the sum which appears in equation (6) by an integral. Using the 
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same techniques developed by Flax and Raich one can obtain the mag- 
netization; 

(s^) = -i- (7) 

' ' 2A 

where 



T 
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K(k) = complete elliptic integral of first kind (13) 

and 



The solution of (8) give as a function of temperature, magnetic field and 
magnetization o Hence by substituting equation (8) into equation (7) the 
magnetization is obtained. 


For an adiabatic process, that is, when there is no change in en- 
tropy, 

•dr = - -L~ dH» (15) 

CRt \ 9t /jjs 


where Cjj is the magnetic specific heat. Since for a ferromagnet 
^9 (s^)/9Tjjjt is negative, an increase in field produces an increase in 
temperature; moreover, the increase will be expected to be largest 
near the critical point. This is known as the magnetocaloric effect. 

The entropy, which is a measure of the order of the system, can be 
calculated from 



The change of magnetization with respect to the reduced temperature 
can be derived from equation (5) as follows: 
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d(s^) ^ 

At jj, 

RESULTS AND DISCUSSION 

In the absence of an applied magnetic field the random phase ap- 
proximation predicts that the magnetization shows a second order tran- 
sition and has a critical temperature which can be written as 



0.3589 


J(0) 


(Ho = 0) 


(19) 


where is the Curie temperature. It is then feasible to relate 
and H’ as follows: 


MggH /j.ggH 

H’ = — X 0. 3589 

j(0) (kgTo) 


(20) 


Figure 1 shows a plot of (s^) as a function of temperature for sev- 
eral values of H for a bcc lattice. When H’ equals zero there are no 
solutions for (s^) above a certain temperature T which is called the 
Curie temperature. This is the point where the long-range order would 
disappear and above this point there is complete disorder of the spins. 
When H' is not equal to zero the ferromagnetic transition occurs not at 
a single temperature, but over a range of temperatures forming a "Curie 
region" in which the transition is smeared out. This is marked by the 
appearance of a "tail" on the magnetization curve. The cause of the tail 
is that the long-range order persists because the spins feel the effect of 
an applied magnetic field. The higher the field the greater the broadening 
of the transition. 
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Figure 2 shows the values of the temperature change in the mag- 
netocaloric effect as a function of temperature. These results were de- 
rived from equations (15) and (18). An expression for Cjj, which is re- 
quired in equation (15) was derived from the expression (0E/0 t)jj, where 
E is the enthalpy. The enthalpy is given by E = U - h(s^). Where U 

is the internal energy of the system, and can be found from a knowledge 

g 

of the Hamiltonian. As indicated in figure 2, the magnetocaloric effect 
is largest in the neighborhood of the Curie temperature. The magnitude 
of the magnetocaloric effect is dependent on the internal magnetic field 
as well as the temperature. For very large fields the magnetocaloric 
effect rises indefinitely as for an ideal paramagnet. (Note that the lattice 
specific heat is not taken into account in this calculation). 

Figure 3 shows a plot of S as a function of temperature for several 
values of H’ . The curves show a characteristic behavior for ferromag- 
netic to paramagnetic transition. It is included here to illustrate that the 
entropy can be calculated for a ferromagnet as a function of temperature, 
magnetization and applied field. Such curves are useful for the design of 
refrigeration cycles. 
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